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Abstract 

i We compute the expectations of the squares of the electric and magnetic fields in the 

' vacuum region outside a half-space filled with a uniform dispersive dielectric. We find 

a positive energy density of the electromagnetic field which diverges at the interface 
despite the inclusion of dispersion in the calculation. We also investigate the mean 
squared fields and the energy density in the vacuum region between two parallel half- 
spaces. Of particular interest is the sign of the energy density. We find that the energy 
^ . density is described by two terms: a negative position independent (Casimir) term, and 

I ' a positive position dependent term with a minimum value at the center of the vacuum 

region. We argue that in some cases, including physically realizable ones, the negative 
04 ' term can dominate in a given region between the two half-spaces, so the overall energy 

^ . density can be negative in this region. 

(N ■ 

^ ; PACS categories: 12.20.Ds, 03.70.-hk, 77.22.Ch, 0A.62.+V. 

o . 

^ ! 1 Introduction 

^: 

I I In 1948 Casimir made the remarkable prediction that there is an attractive force between a 

■ pair uncharged parallel plane perfect conductors Furthermore, he argued that this force 

^ , arises solely from a shift in the energy of the vacuum state of the quantized electromagnetic 

qh| field. An early attempt by Sparnaay [2] to observe this force was inconclusive, but in 

^ • recent years several new experiments 1^, i, ^, |^, Q have been performed which seem to give 

• , good agreement with Casimir's prediction. (To be more precise, most of these experiments 

r> I actually measure the force between a plate and a sphere and incorporate a theoretical 

' correction to compare to Casimir's result. Of the recent experiments, only that of Bressi et 

al uses two parallel plates.) 

If the energy of the vacuum state is zero in the limit of infinite plate separation, then the 
attractive force found by Casimir would seem to imply a negative vacuum energy at finite 
separation. In fact. Brown and Maclay showed that for perfectly conducting plates, 
one has a constant negative vacuum energy density. This conclusion is of great theoretical 
interest, because negative energy density has the potential to cause some rather bizarre 
effects in gravity theory. (See, for example, Ref. Q and references therein.) However, 
questions have been raised as to whether the negative energy density will still arise in a more 



realistic treatment in which the plates are not perfect conductors [10, 11|. In particular. 
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Heifer and Lang [^] calculated the energy density outside of a single half-space filled with 
a nondispersive dielectric material and obtained a positive result. They interpreted this as 
a positive self-energy density associated with a single plate which would add to the negative 
interaction energy density between a pair of plates. Heifer and Lang conjecture that the 
net Casimir energy density might be positive when the self energy is accounted for. If this 
conjecture is correct, then the situation would be analogous to that of the energy density 
in classical electrostatics. A pair of oppositely charged particles have a negative interaction 
energy, but the net energy density, which is proportional to the square of the electric field, 
is always positive. 

However, the Heifer and Lang calculation does not include dispersion, which is essential 



in a realistic treatment. Numerous authors, beginning with Lifshitz |12], have studied the 
effects of dispersion upon Casimir forces. However, these authors have been concerned with 
the force or the total energy, and not the local energy density. The purpose of this paper 
is to present a calculation of the Casimir energy density in a model in which dispersion 
is included. For this purpose, we will use the methods of source theory developed by 



Schwinger and coworkers [13, This is a method based upon the calculation of Green's 



functions which is especially well suited to dissipative materials, and was used by Schwinger 



et al to rederive the results of Lifshitz. Milonni and Shih [15| have used conventional 
quantum electrodynamics to reproduce some of the results of source theory. There has also 
been considerable interest in recent years in quantization of the electromagnetic field inside 
dissipative materials using operator methods The relation between the results 

of the latter set of authors and those of Schwinger et al has not yet been clarified. 

The outline of this paper is as follows: In Sect. ^ we review the source theory approach 
as applied to parallel interfaces of dielectric media. In Sect. ^ we compute the expectation 
values of the squares of the electric and magnetic fields in the vacuum region outside a 
half-space filled with a uniform dispersive dielectric. We extend this calculation to the 
case of two parallel dielectric half-spaces and also discuss the energy density in Sect. 
Conclusions are given in Sect. |^. 

2 Green-Function Approach for Multilayer Dielectrics 

This section is a review of the formalism of Schwinger et al One begins by writing 

the Maxwell equations for the macroscopic electromagnetic fields produced by an external 
polarization source P, which formally describes the zero point fluctuations of the field^ 

B = V X A, 
E = -A - V0, 
VxB = eE + P, (1) 
V • (eE + P) = 0, 

where e is the dielectric constant of the medium. The wave equation for the electric field 
resulting from the Maxwell equations is 

Heaviside-Lorentz units with c = h — 1 will be used in this paper. Also, it is assumed that the magnetic 
permeability is unity. 
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-V X (V X E) - eE = P, 



(2) 



By assuming a linear relation between sources and fields, the electric field can be written 
as a spacetime integral 

B{x) = j d^x'Y{x,x')P{x'), (3) 

where x = {t, r), x' = {t', r') and ^ is a Green's dyadic, which satisfies (^) with a (5-function 
source. Let 



r (r,r',w) = / dre"^^ T {x,x'), (4) 

J — oo 

where r = t—t'. From (|2|) and (^), it follows that ^(r, r', u>) satisfies the following equation: 

-\/ X {V xY)+Lo^eY = -io'^TSir -r'). (5) 



So far, the discussion has been purely classical. At this point, Schwinger et al [|13[| use 
source theory to identify the Green's dyadic F with an "effective product of electric fields" 



'-{E,{T)Ek{v'))=V,k{v^v',Lo). (6) 

We can interpret this as the Fourier transform of the electric field correlation function. 
From the Maxwell equation V x E = — B, one finds the corresponding expression for the 
magnetic field: 



^{Bj{v)Bk{r')) = ej7mefc„p(ViV„//a;^)Fmp(r,r',a;). (7) 

Note that % makes its first appearance in these expressions. These expressions can be 
identified with the vacuum expectation values of products of field operators, which appear 
in the more conventional field theory approach to quantization of the electromagnetic field. 
From now onward, we revert to units in which ^ = 1. In order to calculate the field 
correlation functions, one needs to find the Green's function F occurring in d). This 
amounts to solving a classical boundary value problem. 

The interfaces between the media are chosen to be perpendicular to the z direction, 
so for now it will only matter that the dielectric constant changes in the z direction only. 
Therefore, it is convenient to introduce a transverse spatial Fourier transform 

"r{v,v',u:) = I dk^ ^e'^^(r-r')^T{z,z',k^,u;), (8) 



3 



where the vector k_|_ can be chosen to point along the +x axis (k = |k_|_|). 



Some components of T are found to be |13] 
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where e' = e(z'), and g^ , the "transverse electric", and g^ , the "transverse magnetic" 
Green's functions satisfy 



did k^ 
dz e dz e 



ijj 



g^'iz^z') 



5{z - z), 
5{z-z'). 



By introducing the quantity 



(|lO|) can be written as: 



2 7 2 2 

K, = k — uj e, 



2 

dz^ 
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(10a) 
(10b) 

(11) 

(12a) 
(12b) 



So, in order to find the field correlation functions as defined in (P) and (Q) in a given 
situation, one needs to solve these equations with the appropriate boundary conditions. We 
consider here two cases. 

3 One Interface Case 

We now specialize the above discussion to a situation in which the inhomogeneity of the 
dielectric constant is due to a plane interface separating a dielectric substance from a 
vacuum: 



z > : e(z) = 1, 
z < : e(z) = e^. 
Here is a function of frequency, but not of position. 



(13) 
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3.1 Boundary Conditions 

In solving ( |12aD and ( |12bD , we use the following boundary conditions. At z = z' , g is 
continuous but the derivative is discontinuous at this point [|l9|]: 



dz 



(14) 



At the boundary (z = 0) we use the conditions for continuity of E^, Ey, eEz, and Bi. The 
first three, as seen from (|6|), imply the continuity of T^x, ^yy, and eV^z and subsequently, 
from (^, the continuity of , g^ , and 

1 d I d s 
€ dz e' dz' ^ 

The continuity of Bx implies that of VzV^'Fj^y , as seen from Eq. (p4|), which is given below. 
From this, using (^) and (pb|), we deduce the continuity of dg^ /dz. 
The solutions g^ and g^ in the vacuum region have the form 



" = 2;^ • <'^''' 

„-Ko\z-z'\ I „lp-KQ{z+z') 

9^ = , (15b) 



where 



1^0 - 



(16a) 
(16b) 



Here kq and ki represent the quantity k as defined in (11) for the vacuum region (e = 1) , 
and for the dielectric half-space region (e = e^) , respectively, and r and r' can be identified 
as reflection coefficients for two polarization states, _L and || respectively, corresponding to 
electric fleld vector being perpendicular or parallel to the plane of incidence of an linearly 



polarized electromagnetic wave [|19 



3.2 The Electric Field 

Using (^), we write the formal expectation value of the square of the electric fleld at coin- 
cident points as 

/oo roo -j^ ^ roo roc 
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In the second step, we assumed that the integrand is an even function of uj. By complex 
rotation (u this becomes: 



oo roo 



;^'>^ = ^ / dC I dkkVkk- (18) 



Jo 



Note from (11) that > when w is imaginary. By means of all of the components 
of r in a given region can be written in terms of Txx and Tyy-. 

Tzx{z,z') = -^-^Txx{z,z'), (19) 

^ , k"^ d d ^ . uj"^ ^, 

T^^{z,z) = -^ — —Txx[z,z) + ^d{z-z). 

(k^J OZ OZ' 

By taking the limit z z', and thus omitting the delta function, T^k becomes 



[K ) 



or by (I), using e = 1, 



kk 



1^2 

+ V,V,.g^ + -^V,V,,(V,V,/5^) 
= wV + (A;' + V,Vy)5^. (21) 
Using ( pA[ ) and (15), this becomes 

rkk = - + :^ [io'r + {2k' - co') r'] e'^-. (22) 



Equation (18) gives a formal expectation value only, because the integral is divergent. 
However, the divergence comes only from the lo'^/k term in and is independent of z. It 
is the usual empty space vacuum divergence. We will henceforth drop this term and denote 
the resulting finite expectation value by {E'^)- The renormalization results in a quantity 
which vanishes at large distances from the interface: {E'^) — > as z — > oo, which amounts 
to finding the difference in (-E^^ with the boundary and without it. Thus we find 



1 poo poo u 

{E') = ^J dCj^ dk-^ [-CV + (2A;2 + C^) r'] e'^- 



(23) 
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3.3 The Magnetic Field 

Now we compute the expectation value of the magnetic field. Using Eq. (^), we find 



i{By{r)By{r')) = ^ (V^V./F^^ - V^VyF,^ - V^V^/F^^ + V^V^/r,,) , (24) 



i(S,(r)5,(r')) 



1 



{^y^y'^xx ^ y^ x'^xy ^ x^ y'^yx ~t~ ^ x^ x'^yy^ 



From the definition of k_L, it follows that all derivatives in y vanish, so we can write the 
sum of the above terms as 



i{Bi{r)Bi{r')) = {V.V.'Tyy + V^V,/F. 



VxV^'F^^. — VzVx'^xz + V^^Vj^/F^^ + Vx>Va;/F, 



Using (19), we have 



(25) 



Va;Vx'F;2^(r,r ,a;) 
V^Vj:'F^^(r,r',u;) 
V^.V,/F,^(r,r',w) 



This leads to 

i{B{x)B{x')) 
Using (P), this becomes 

i{B{x)B{x')) 



f dio 


f dk± 


" 1 




(2vr)2 





dkx (A;' 
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V,V,'F, 



(2^)2 (k2)2 

dkx A;2 



(27r)2 AC2 



(27r)2 k2 



^z^^'Fxa:, 
^ z''^ XX • 



(/c2 + V,V,,) Fj^j,(z, z') + _V,V,/F^^(z, z') 



/ ^ /^[(^^ + v.v.,).^(.,.')+-V(.,.')] 



(26) 



(27) 



(28) 



Following the same procedure as used above in calculating {E"^), we find the finite mean 
squared magnetic field to be 



(29) 



Note from ( p3|) and ( |29| ) that {E"-) <^ {B'^) under interchange of r and r'. Now, the mean 
energy density can be calculated as 
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(30) 



Using (|23D and (|29|), this becomes 



1 /-oo foo 7,3 



(31) 



We can write ?7 in a form more convenient for numerical calculation by introducing polar 
coordinates u and ^ (C = ucos9, k = usinO): 



U 



47r2 







duu^ / de{sm9f (r + /) e 



-2uz 







We use the Drude model for the dielectric function 



(32) 



1 



where cOp is the plasma frequency. From (16a), ( |16b| ), and (p^), we find 



(33) 



By a substitution (cos(0) 



u — A / li^ + 



u^{cos9y +ujp- uicosOyJv? + lo'^ 



^^(cos 9)'^ + Up + u(cos 0)^^/ u'^ + 
t), U becomes 



(34a) 
(34b) 



U = —J duu^ J dt {I - t"^) {r + r') 
By the same coordinate transform, (|2^ ) and (^) become 



-2uz 



(35) 



1 

4^ 
1 

4^ 



-2uz 



-2uz 



(36a) 
(36b) 



The plot for (-E^) and (i^^) , as well as U is shown in Figure 1. As we can see from the 
figure, the energy density is positive. Now we consider some limiting cases. 
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Figure 1: 

The expectations of the squares of electric field, magnetic field, and energy density near 
the dielectric half-space are illustrated. 

3.4 The Fields near the Interface 

To see how U behaves for small z (large u) , we first perform the t integration in ( |35|) , which 
can be done analytically, and then Taylor expand the resulting expression in the brackets 
in powers of u~^. That is, we are expanding all of the integrand except for the exponential 
factor. To the leading order we find: 



The asymptotic behavior of the mean squared fields ( |36a| ) , ( ^6b ) in this limit is 

(^=> ^ #i (38a, 

We see that {E'^) dominates over (B^), so that U 1/2 {E"^); this is due fact that the 
leading order in the expression in braces in ( p6a| ) is oc as compared to u~'^ in ( p6b ). If 



we compare these expressions to ones that would result if dispersion were not included in 
the calculation, it can be seen from ( |36| ) that in this case {E^'j oc z~^, and same for (^B"^ 



(see also [|^). As seen from (38), the inclusion of dispersion in the calculation reduces the 
power in z up to two orders, but it does not remove the singularity of the results at z = 0, 
as might be naively expected. 

After more careful consideration, it is not surprising that dispersion alone is insufficient 
to render the results finite at the interface. The integrals for (i?^) and (-B^) at z = 
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diverge quartically in a frequency cutoff. However, in general dielectric functions go as 



e{uj) ~ 1 + 0{lo-^) (38c) 

as a; ^ oo. Thus the reflection coefficients will go to zero no faster than a;~^, leaving the 
integrals quadratically divergent. This argument explains why (-B^) oc z~'^ for small z, but 
understanding the behavior of (-E/^) requires examining the dependence of the reflection 
coefficients r and r' upon the transverse momentum k. In fact the contribution of the coef- 
ficient r, which describes modes with the polarization vector perpendicular to the plane of 
incidence, does go as z~'^. This coefficient depends only upon frequency, and falls as w"^ for 
large uj, as can be seen from ( |34a| ) . The coefficient r' describes modes with the polarization 
vector parallel to the plane of incidence, and goes to one as 9 7r/2 (corresponding to 
grazing incidence) for all frequencies. It is this behavior which leads to the z~^ singularity 
in (-E^) and hence in U. (The role of cutoffs for the quantized electromagnetic field in 
dielectrics has been discussed in more detail by Candelas ||2^. Barton ||2l| has recently 
emphasized the fact that dispersion alone will not remove all divergences.) 

The divergence of U is not considered to be physical, but as resulting from the ideal- 
ization of the wall as a perfectly smooth surface. One way of removing this singularity is 



to allow the position of the boundary to fluctuate |22]. It seems plausible that such effects 
as surface roughness, or the atomic nature of matter on small scales can also introduce a 
physical cutoff that makes the mean squared fields and the energy density finite everywhere. 

3.5 Case of a Perfect Conductor 



Now we consider the limit e ^ oo. In this limit, as seen from (|TT|), ( 16aD and (fBB), r ^ — 1, 



and r' —>■ 1 . Equation (|3ll) implies that U becomes zero, as expected, and (|36aD and ( |36bD 
give: 

(B^) ~ -I^.^. (40) 



These well-known results are consistent with the asymptotic Casimir- Polder potential P3| : 



where ao is the static polarizability of an atom near the interface. 

4 The energy density between dielectrics 

In this section we calculate the energy density in a vacuum region of width a between two 
dielectric half-spaces. We define the dielectric constant as: 
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0<2;<a:e = 1, 
z < and z > a : e = e^. 



(42) 



In the vacuum region, occurring in ( |12a| ) has the form 

^ If exp(— kq I z — z' I) + exp(— 2Koa) exp(Ko I — -z' |) 



+ 



2kq [ exp(— 2Koa) — 1 

r [exp(— Ko(z + z')) + exp(— 2Koa) exp(K(z + z'))] 
exp(— 2Koa) — 1 



+ 



(43) 



and has the same form as g^' with r and r' interchanged, where r and r' are defined in 
(|l6al) and ([l6b|) . We can now calculate the electric and magnetic fields in this region. 



4.1 The Electric Field 

The first and the second term on the right hand side of (^) for the present case, in the 
limit z ^ z' , are 



SO that r^fc becomes 



a;2 1 + ^2g-2Ka 

2k 1 — r^e^^'^" 

1 + r'^e" 
2k 1 - r'2e-2«« 



2 ^ f' ^ 2 j ^ 2 Aijtx ^2 Aij/S 
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+ K^ /(e-2K2 _)_ e^""""e' 



2k 
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(44) 
(45) 



fcfe 
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J2 



+ 



g2/^(i ^^/2 



+ (2A;2 _ c^2) 



+ 



1 — r 



l1n—2Ka 



e-'^'^ cosh [k(2z - a)] . (46) 



We again drop the first term on the right hand side of the above expression. After intro- 
ducing polar coordinates and using ( [ll| ) and (^6|), ( |T8|) we find 



1 



duv:"- I de sin((9) <^ cos2(6') 



+ 



cos2(^) 



J2 



+ 



+ 



1 — r^e 



2„— 2«a 



+ (1 + sin- 



j,/2 g—2ua 



■ cosh [u {2z -a)]} , (47) 



or with (cos(0) — > t) 



,/2 



+ 



+ 



+ 



1 — r^e 



2n—2ua 



+ 2 



j,/2g— 2ua 



e""'' cosh [-u (2z - a)] 



(48) 
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Figure 2: 

The expectation values of the squared electric and magnetic fields, as well as the energy 
density in the vacuum region, are illustrated for Upa = 200. 



4.2 The Magnetic Field 

In the same way, (|2^) leads to 



1 

2^ 



+ 



Jo 

r 



^'2 



+ 



+ 



1 — r'^e 



l2„-2ua 



e-'"" cosh [u {2z - a)] 



(49) 



This expression differs from the one for (-E^) only in the z- dependent term with r ^ r' . A 
plot of (-E^) , (^^) and U is shown in Figure 2 for uopa = 200. It can be seen in this figure 
that significant, but not complete cancellation occurs between (-E^) and (5^) . 



4.3 The Energy Density 

Now using (^0|), (^), and (^9|), the energy density in the vacuum region can be calculated 
as 



U 



1 

2^ 



duu 



dt < t 



^/2 



c2ua 



+ 



,/2 



o2Ma 



+ 



1 



+ 



1 



j./2g— 2ua 



e-^'^ cosh [n (2z - a)] 



(50) 



By analyzing this expression we can make some conclusions about the sign of U . First we 
note that it is position dependent, and we also note that the first term is always negative 
and the second term is always positive. The overall sign of U depends on the choice of a 
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0.3 




Figure 3: 

The energy density in the vacuum region between two dielectric half-spaces is illustrated 
for three values of the parameter iVpa. The dashed horizontal line is the energy density for 
the perfectly conducting limit, (|5^). 



and ujp. As ujpU grows, U at the midpoint decreases, becoming negative for Upa 100, as 
seen in Figure 3 and Figure 4. 

In Figure 4 we see how the energy density at the center of the vacuum region changes as 
the product copa increases. It can be seen both in Figure 3 and Figure 4 that U approaches 



the value given in (52) as ujpa becomes large. The separation at which U becomes negative 
at the center is 

99 /14.8eV\ 
a > Qc = — = l.o/^m , (51) 



where 14.8eV is the plasma frequency of aluminum. 



4.4 A Perfect Conductor Case 

In the limit ^ oo, r — > — 1 and r' — > 1. Then only the first (z-independent) term survives 
in (|50[), and we get the familiar result |^: 



U 



■K 



' 720a4 ■ 



(52) 



It is also interesting to examine the expressions for (E"^) and (B"^) in this limit. After 



performing the t integration, (48) can be written as 



1 

2^ 



du u 



3— 2u(a— z) 



3 1 



o2«a 



+ 



1 



-2ua 



+ 



-2ua 



(53) 
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Figure 4: 

The graph represents the energy density at the center of the gap between the two dielectric 
half-spaces as a function of ujpa. As seen from the graph, the energy density at center 
becomes negative when copa ~ 99 or larger. Again the dashed line is the perfectly conducting 
limit. 

and after performing the u integration as 



vr^ 1 

+ 



^(3)ff)+^(3)Ci_f 
Va/ \ a 



(54) 



720a4 327r2a4 

Here ip^^"^ = -4^ lnr(2;) is the polygamma function of order three. It satisfies the reflection 



formula ||2J 



^(3) f£)+^{3)M_~ 



( z\ 



This yields 



vr 



+ 



vr^ l + 2cos^(7r 



In the same way, one finds 



720a4 16a4 sin^'fvr^l 



7r2 l + 2cos^(7rf) 



720a4 16a4 sin^(7rf) 



These results are in agreement with those given by previous authors p5|, 2t\ 



(55) 



(56) 



(57) 



4.5 Energy density near the Boundary 

To see how U grows near the interface we note that for small z (large u), (|5^) becomes 
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1 1'°^ 

J duu^ J dt {1- t^) (r + /) e-2«^, (58) 

so it reduces to (|35[), the solution for one interface case. In the z — > Umit, this expression 
reduces to (^7|) . By the same reasoning, the expressions for (-E^) and {B"^), (|48|) and (^9|) 
reduce to (|36a| ) and ( |36bD , respectively, in the small z limit. 



5 Conclusion 

In this paper we have examined the mean squared fields and the energy density in the 
region between a pair of half-spaces filled with dispersive media. We found that these 
quantities diverge at the boundaries of the media, despite the inclusion of dispersion in the 
calculation. This divergence indicates a breakdown of a continuum description in which 
the dielectric function changes discontinuously at the boundary. It also shows that there is 
a positive self-energy density in the region outside of a single plate. Nonetheless, we have 
also found that it is possible for the net energy density in the region between the plates to 
be negative, depending upon the plate separation and the plasma frequency of the material 
involved. The existence of an attractive Casimir force is not an indicator of whether the 
energy density at the center of the plates is actually negative or not. 

We have found that the energy density at the center becomes negative when ujpa > 
100. Thus for fixed plasma frequency LOp, the energy density always becomes negative for 
sufficiently large separation a. Of course, in this limit the magnitude of the energy density 
is also becoming small. Similarly, for fixed a, the energy density becomes negative for 
sufficiently large LOp. In the limit that uj p ^ oo our results approach the constant negative 
energy density of the perfectly conducting plates. It should not come as a surprise that 
there is a regime of negative energy density. The calculation assuming perfect conductivity 
does have a region of validity so long as ujp is large and one is not too close to one plate. 
Qualitatively similar behavior has recently been found for the vacuum energy density near 
a domain wall p7| . 

In this paper, we assumed a particular form for the dielectric function, (p3[), given by the 
collisionless Drude model. This is a good model for many metals especially alkali metals, 
and is the generic form for the dielectric function of all materials at high frequencies. Thus 
taking a different form for Cd would change the details of our results, especially far from 
an interface, but should lead to the same limiting forms near an interface. We have also 
assumed zero temperature throughout this paper. For systems at room temperature, this 
should be a good approximation when the separations are of the order of a few micrometers. 
More generally, one can ignore thermal effects at distances small compared to \/{kT). 



Contrary to the view expressed by Lamoreaux |11], the appearance of negative energy 
density in a quantum field theory is very natural. One can easily find quantum states of 
the free quantized electromagnetic field in empty space which have local negative energy 



densities. A squeezed vacuum state is an example |28, 29 1. The energy density of a quantized 
field has to be defined as a difference between that in empty Minkowski spacetime, and that 
in a given state and is no longer positive definite, as it was for a classical field. Apart from 
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coupling to gravity, which produces extremely small effects, no clear way has been found 
to directly observe the local energy density. In certain limits, the negative energy density 
in a squeezed vacuum state has been shown theoretically [29| to produce an effect on the 
magnetic moment of a spin system. Whether this effect could ever be observed, and whether 
negative Casimir energy density can produce similar effects is unknown. 
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